STOCHASTIC DOMINATION AND COMB 
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ALEXANDER E. HOLROYD AND JAMES MARTIN 

Abstract. There exists a Lipschitz embedding of a d-dimensional 
comb graph (consisting of infinitely many parallel copies of Z d_1 
joined by a perpendicular copy) into the open set of site percolation 
on Z d , whenever the parameter p is close enough to 1 or the Lips- 
chitz constant is sufficiently large. This is proved using several new 
results and techniques involving stochastic domination, in contexts 
that include a process of independent overlapping intervals on Z, 
and first-passage percolation on general graphs. 



1. Introduction 

The following natural generalization of percolation theory is prompt- 
ed by the results of |2], [/J. Let G and H be graphs. For p e [0, 1], 
consider the site percolation model on H, in which each vertex is open 
with probability p, and otherwise closed, independently for different 
vertices. An embedding of G in the open set of H is an injective 
map from the vertex set of G to the set of open vertices of H, such 
that neighbours in G map to neighbours in H. Define the critical 
probability 

p c (G,H) := 

inf|p : P(3 an embedding of G in the open set of if) >0\. 

If Z + is a singly-infinite path then p c (Z + ,H) is simply the usual 
critical probability p c (H) of site percolation on H (see e.g. |5j for back- 



ground). For the doubly-infinite path Z, it was proved in [11, Proof 
of Theorem 3.9] that p c (Z, H) also equals p c {H) for any infinite con- 
nected H. Observe that if G, H are subgraphs of G', H' respectively 
then Pc (G,H')< Pc (G',H). 
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Figure 1 . Part of the comb graph K in dimension d = 
2 (left) and d = 3 (right). 

We focus on the question: for which graphs is it the case that 
p c {G, H) < 1? Let Z d be the usual cubic lattice, with vertex set 
also denoted Z d , and with vertices x, y joined by an edge whenever 
\\ x ~ y\\i = 1- Also let Zm^-i denote the spread-out lattice, in which 
vertices x,y G Z d are joined whenever < \\x — y\\oo < M. It was 
proved in and |7| respectively that p c (Z d_1 ,Z^) < 1, while on the 
other hand p c (Z d , Zf M] ) = 1 for all M. 

An embedding of G into Zj^j may also be regarded as an M-Lipschitz 
embedding of G into Z d . In that language, the results mentioned in 
the previous paragraph say that M-Lipschitz embeddings of Z d_1 into 
Z d are possible whenever p or M is large enough, while Lipschitz em- 
beddings of Z d into Z d are never possible for p < 1. 

In this article we address lying between the last two mentioned 

above. Define the d- dimensional comb graph K d to have vertex set 
Z d , and edges (z, z + e$) for every z G Z d and all i with 1 < i < d — 1, 
together with (2, z + e^) for all z such that z% — (where ei, . . . , 
are the standard basis vectors). Thus, K d consists of a stack of parallel 
copies of Z d_1 (perpendicular to coordinate d), connected by a single 
perpendicular copy of Z d ~ x (perpendicular to coordinate 1). For d > 2, 
K d is isomorphic to the product of the 2-dimensional comb IK 2 with 
Z d ~ 2 . See Figure [l] for illustrations of K 2 and K 3 . 

Theorem 1 (Comb percolation). We have p c (K d , Z^) < 1 for all 
d>2. 

Corollary 2. For all d > 2 we have p c {K d , Z d M] ) as M ->■ 00. 

Our proof gives an explicit upper bound for p c (K d , Z^i), but we have 
not attempted to optimize it. The spread-out lattice ZjL in Theorem [I| 
cannot be replaced with the nearest-neighbour lattice 7L d . Indeed, it 
was proved in JT] that j9 c (Z 2 ,Z d ) = 1 for all d > 2; since Z 2 is a 
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subgraph of K. d this implies p c (K d , Z d ) = 1 for d > 3. It is also easy 
to see that p c (K 2 , Z 2 ) = 1, since the backbone of K 2 would have to be 
embedded as a straight line in Z 2 . On the other hand, our techniques 
may be adapted to prove p c (K , H) < 1 for some graphs H with edge 
sets intermediate between those of Z d and ZjL - in particular it seems 
plausible that this could be done for the "star lattice" Zj^, but we 
have not pursued this. Such questions reflect details of the local lattice 
geometry, whereas the fact that p c {K d ,Z d M ^) < 1 for large enough M 
(as implied by Theorem [TJ) is more fundamental. 

Our proof of Theorem [l] will make use of several new results and 
techniques involving stochastic domination, which we believe are of 
independent interest and wider applicability. Stochastic domination 
by i.i.d. processes is a powerful technique for proving results of this 
kind, because it enables facts proved for the i.i.d. case to be transferred 



to other settings. One widely used tool is the result of JT0[ that a k- 
dependent Bernoulli process with sufficiently high marginals dominates 
any given i.i.d. product measure. However, the key process that we will 
need to control (of "bad points") is not /c-dependent, and in fact is not 
dominated by any product measure. Therefore the methods we use are 
of a different nature. 

Background on stochastic domination may be found in ||, Ch. II, §2], 
for example. For our purposes, the following definition via coupling 
will suffice. Let X and Y be random variables taking values in the 
same partially ordered space. Then we say that X stochastically 
dominates Y if there exist X', Y' on some probability space with X' 
and X equal in law, Y' and Y equal in law, and X' > Y' almost surely. 
The underlying partial order will be inclusion (in the case of random 
sets) or pointwise ordering (for real functions). 

Our first tool is a simple but useful stochastic domination result on 
overlapping intervals in a one-dimensional setting. For c G (0, 1), say 
that a random variable X has geometric distribution with param- 
eter c, denoted Geom(c), if F(X = r) = (1 — c)c r for r = 0, 1, 2, ... . 
(Note that the value is included, and that c is the probability of a 
"failure" rather than a "success"). In the following, the interval (a, b) 
is taken to be empty if a = b. 

Theorem 3 (One-dimensional domination). Let (G n ) ne z be i.i.d. 
Geom(c) random variables. The random set Z fl U n( z Z (n — G n , n + G n ) 
is stochastically dominated by the open set of i.i.d. site percolation on 
Z with parameter min(A^c, 1). 

Our second tool concerns first-passage percolation. As we explain 
in Section 0, it can be regarded as unifying and generalizing ideas in 
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|], U, ||. Let V be a countable vertex set. For every pair of distinct 
vertices x,y G V, the directed edge e = (x,y) is assigned a random 
passage time W(e) = W(x,y) taking values in [0,oo]. The passage 
times of different edges are independent but not necessarily identi- 
cally distributed. (We can model a process on a graph other than 
the complete graph by taking some passage times to be oo almost 
surely.) In addition, each vertex x G V has a deterministic source 
time t(x) G (—00,00] at which it is "switched on". (For example, to 
model growth started at a single source a we would take t(a) = and 
t(x) = 00 for all other x.) The occupation time of x G V is the time 
it is first reached: 

Vm = X < k = l ' 

We now consider a collection of countably many models on the same 
vertex set, indexed by % G X. Different models have identically dis- 
tributed passage times, and are independent of each other, but may 
have different source times. Write U(x) for the source time of x in 
model i, and Ti(x) for the occupation time of vertex x in this model. 
Let 

T(x) := inf Ti(x). 

i 

Finally, consider another model with source times given by 

t(x) := inf ti(x), x G V, 

i 

and with the same passage time distributions as the other models. 
Write T(x) for the occupation time of x in this model. 

Theorem 4 (First-passage percolation domination). Under the above 
assumptions, (T(x)) xe v ^ s stochastically dominated by (T(x)) xe y. 

We prove Theorems |^ and § at the end of the article. In the next 
section we explain how these results are used in the proof of Theorem |l|. 

2. Outline of Proof 

In this section we explain the main ideas behind the proof of Theo- 
rem |1|. Our starting point is the following strengthening of a result of 
H (the latter has been applied in pi ^J, and extended in other direc- 
tions in M). For x — (xi, . . . , x^-i) G Z d_1 and z G Z we denote their 
concatenation thus: (x,z) := (xi, . . . , Xd-i, z) G Z d . Vertices of Z d will 
sometimes be called sites. 
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closed sites; 

stacked Lips- 
chitz surfaces L n 
avoiding closed 
sites; 

good sites 
(where L n is as 
low as possible); 

perpendicular 
Lipschitz surface 
H avoiding bad 
sites; 

three selected 
sites y; 

obstacles A y at 
those sites y. 



Figure 2. The main objects used in the construction of 
the embedding. Here d = 2 and p = 0.885. 



Theorem 5 (Stacked Lipschitz surfaces). Consider site percolation on 
Z d with d > 2. If the parameter p is sufficiently close to 1 then a.s. 
there exist (random) functions L n : Z d_1 — >■ Z, indexed by n G Z ; tuii/i 
i/ie following properties. 



(la) 
(lb) 

(lc) 
(Id) 



The site (x, L n (x)) G 7L d is open for all x G Z d_1 and n G Z. 
For eac/i n, i/ie function L n is 1-Lipschitz in the sense that 

\L n {x) — L n (x')\ < 1 whenever \\x — x'Hoo = 1. 
L n (x) > 2n for all x and n. 
L n _i(x) < L n (x) for all x and n. 



For each n, the graph L n (x)) : x G Z d_1 } of L n is a "Lipschitz 
surface" , and Theorem [5] asserts the existence of an ordered stack of 
disjoint open Lipschitz surfaces. See Figure 0. This strengthens the 
result of that one such surface exists for p sufficiently close to 1. Our 
Lipschitz surfaces differ from those in 0, || in that we use the oo-norm 
rather than the 1-norm in (|T5]) - this is relatively unimportant, but will 
be convenient for our construction. The condition (flc| ) will be helpful 
in keeping track of the typical position of each surface. 
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In proving Theorem |5| we will define a particular family of functions 
(L n ) having additional desirable properties. In fact, (L n ) will be the 
minimal family satisfying (p~a|)— (pT5|) in the sense that for any other such 
family (L' n ) we have L n (x) < L' n {x) for all x,n. 

Our aim is to weave these Lipschitz surfaces together using another 
Lipschitz surface perpendicular to the stack. Observe that the mini- 
mum possible value of L n (x) is 2n + 1. We pay particular attention 
to those positions where this minimum is attained. Let Z even be the 
set of even integers and Z odd the set of odd integers. We call the site 
(x, 2n + 1) G Z d_1 x Z odd good if L n (x) = 2n + 1, and otherwise bad. 
Note that these definitions apply only to sites whose last coordinate is 
odd, and depend on the choice of the functions L n . Since (x, L n (x)) is 
always open, every good site is open. 

Theorem 6 (Perpendicular Lipschitz surface). Fix d>2. For p suf- 
ficiently close to 1, the functions L n of Theorem [| may be chosen so 
that almost surely there exists a function H : Z d ~ 2 x Z odd — > Z with the 
following properties. 

(2a) The site (H(u),u) is good for all u G 7L d ~ 2 x Z odd . 
(2b) \F[(u) — H{u')\ < 1 whenever: 

\ua-i — u' d _ x \ < 2, and \ui — u[\ < 1 for i < d — 2. 

See Figure 0. The set {(H(u),u) : u G Z d ~ 2 x Z odd } forms a kind 
of Lipschitz surface perpendicular to the 1 coordinate direction. (The 
"<2" in ( pb| ) reflects the appearance of Z odd in the domain of H. Note 
that the d — 1 coordinate of u becomes the d coordinate of (H(u),u).) 
It is relatively straightforward to check that any functions L n and H 
satisfying (pi|)-(|Tcl|) and (|2a|)— (pLj|) give rise to an embedding of K d in 
the open set of ZfL, as required for Theorem |T[ This is verified in 
Section the function H gives the backbone of the comb, while the 
L n 's give the fins. Therefore our main task is to prove Theorems || and 

! 

We will prove Theorem |^ via an extension of the methods of || : the 
Lipschitz surfaces will be constructed as duals to paths of a certain 
type, called A-paths. Now, since the property (pbD required for H is 
essentially property (|TT5|) of our Lipschitz function L (modulo a change 
of coordinate system), an appealing idea is to try to deduce Theorem | 
from Theorem [|. The problem, of course, is that the process of good 
sites is not i.i.d. It is also not dominated by any i.i.d process (because a 
vertical column of k consecutive closed sites gives rise to a bad set with 
volume of order k d ). Nonetheless, we will indeed deduce Theorem ^ 
from Theorem [|, using stochastic domination in more subtle ways. 
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We will proceed by re-expressing the process of bad sites. For each 
y G Z d_1 x Z cvcn we will define a random finite set A y , called the 
obstacle at y, in such a way that 



The field of obstacles will have the stationarity property that (A y + z) y 
is equal in law to (A y+Z ) y for all z. The obstacle at y will be the 
set of points that can be reached from y by A-paths satisfying certain 
conditions. 

The random sets A y will not be independent of each other (since 
the paths used in their construction are shared between different y's). 
However, we will prove that they can be replaced with independent 
sets in the following sense. Let {A y ) y be mutually independent random 
sets, with A y equal in law to A y for each y. We will show 

(3) \ \ A y is stochastically dominated by I \A y . 



A similar fact was proved in || in the context of a Lipschitz percolation 
model. An analogous property for a continuum percolation model was 
obtained in [|IJ , and related ideas appeared earlier in [|J . We will prove 
(|3]) by expressing A y in terms of a first-passage percolation model (via 
the paths involved in its definition), and appealing to the much more 
general Theorem |j. 

Our task is now reduced to proving the existence of a Lipschitz sur- 
face H (as in (|2b~|)) avoiding a collection of independent sets A y . The 
ideas behind the proof of Theorem [| will easily show that the radius 
around y of the random obstacle A y (and thus A y ) has exponential tails 
for p sufficiently close to 1. However, for d > 2, this is not enough to 
allow domination of A/ by an i-i-d. percolation process, since there 
the probability of a closed ball of radius r decays exponentially in r d . 

The final ingredient is a deterministic observation which allows us 
to reduce to a one-dimensional process and hence overcome the above 
dimensionality problem. Here it is important that the object we seek 
is a Lipschitz surface. For x G Z d and r > define the ball B(x,r) := 
{z G Z d : \\x — z\\oo < r} and the one-dimensional stick S(x, r) := 
{x + aed : a G Z and \a\ < r}. 

Lemma 7 (Balls and sticks). Suppose h : Z d_1 — > Z is 1-Lipschitz (i.e. 
\h(x) — h(x')\ < 1 whenever \\x — x'\\oo < 1)- If the graph {(x,h(x)) : 
x G Z^ 1 } does not intersect the stick S(y,2r — 1) then it does not 
intersect the ball B{y,r). 



{x : x is bad} 




y 



y 



s 
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Figure 3. A Lipschitz surface avoids a ball provided it 
avoids a stick. 

See Figure |3] for an illustration. Using Lemma [7], it suffices to con- 
struct a Lipschitz surface that avoids a union of sticks [J xeIi d S(x, G x ) 
with i.i.d geometric sizes G x . This union consists of independent one- 
dimensional processes in each vertical line. Therefore we can use The- 
orem [5] to dominate it by an i.i.d. percolation process (with parameter 
that tends to as p — » 1), and deduce Theorem [| from Theorem |5|, 
and hence complete the proof of Theorem [I]. 

In the next four sections we carry out the steps outlined above to 
prove Theorem [TJ. The stacked surfaces L n are constructed in Sec- 
tion [| Obstacles are defined and dominated by independent sets in 
Section f|, and their radii are bounded in Section [|. The remaining 
details (including the stick argument) are completed in Section |. Fi- 
nally we prove the general domination results, Theorems § and f|, in 
Sections |8| and [7| respectively. The first-passage percolation result is 
proved via dynamic coupling. For the one- dimensional domination re- 
sult we employ a queueing interpretation. 



3. Stacked Lipschitz surfaces 

In this section we prove Theorem [|. We first construct the functions 
L n , and then prove that they have the required properties. We some- 
times refer to the positive and negative senses of the d coordinate as 
up and down respectively, and the other coordinates as horizonal. 

Define a A-path to be a sequence of sites z(0), z(l), . . . , z{m) G Z d 
such that for each i < m, 

(4) z(i + 1) - z(i) G {e d } U A, 

where 

A := I -e d + ^ a i e i : («i, « 2 , • • • , «d-i) G {-1, 0, 1} 
^ i=i 
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That is, each step is up or down, but the down-steps may also be 
diagonal; there are 3 d-1 different types of down-step since each of the 
first d—1 coordinates is allowed to remain the same or change by 1 in 
either direction. (Our definition of a A-path differs slightly from that in 
0, where only 2d + 1 types of down-step were allowed. The difference 
reflects our use of the oo-norm in (|lb|).) For an integer r > 0, we call 
a A-path r-open if its up-steps have distinct locations, and at most r 
of them end with an open site, i.e. among the indices i < m for which 
z(i + 1) — z(i) = td-, the sites z(i + 1) are all distinct, and at most r of 
them are open. We write y — > z if there is an r-open A-path from y to 
z. 

Now define the random set of sites S n by 

(5) S n := {z : y A z for some r and some y with = 2(n — r)}. 
Then let L n be the function whose graph lies just above S n : 

(6) L n (x) := min{£ G Z : (x,£) G" S n }, 
(where min0 := oo). 

Proposition 8. Let the functions L n be defined as above. If p is suf- 
ficiently close to 1, then a.s. L n (x) < oo for all n and x, and the 
properties flla]) -(|ld|) in Theorem [| all hold. 



Proof. From the definition and the underlying stationarity of the per- 
colation process, the process (L n (x)) njX is stationary in the sense that 
(L n+ k{x + y) — 2k) n ^ x has the same law for any k G Z and y G Z d_1 . 
Hence for the first claim it is enough to show that -Lq(O) < 00 a - s - 

In fact we will show that L (0) has exponential tails. For h > 0, we 
have P(L (0) > h) = P((0, h) G So), and this is at most the expected 
number of r-open A-paths from the hyperplane Z d_1 x {— 2r} to (0, h), 
summed over all r. For such a path, let C be the number of up-steps 
that end in a closed site, let U be the number of up-steps that end in an 
open site, and let D be the number of down-steps (including diagonal 
steps). Since the path is from Z d_1 x {— 2r} to (0,h) we must have 
C + U - D = h - (-2r), i.e. U = D- C + h + 2r. Since the path is 
r-open we have U < r, or equivalently A > where A :— r — U. 

For given U,C,D, the number of ways to choose a A-path ending 
at (0, h) together with an assignment of states open and closed to its 
up-steps is at most K u+C+D , where K := 3 d_1 + 2. (There are 3 d_1 
possible directions for a down-step, and two possible states for an up- 
step). For any such choice, the probability that the chosen states match 
the percolation configuration is p u q c < q c , where q :— 1 — p. 
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Therefore 



P(io(0) >h)< K u+C+D q c 



U,C,D,r>0: 
C+U-D=h+2r, 
U>r 



^D+h+r+A 



< J2 K 2D+h+2r q l 
A,D,r>0 



A>0 D>0 r>0 

which converges (exponentially fast) to as h — > oo whenever q < K~ 2 . 
(For the second inequality above, we rewrote U and C in terms of A 
and dropped the conditions U > and C > 0.) 

Now we verify properties ([la|)-(]Id]). For (pTa|) , observe that, for 
some r as in the definition of S n , there is an r-open path to the 
site (x,L n (x) — 1), but there is none to the site (x,L n (x)). Thus the 
site (x,L n (x)) must be open - if it were closed, the r-open path to 
(x,L n (x) — 1) could be extended one step upward (or else it already 
passed through that site). 

Next, note that from the definition of S n , if z G S n then z + v G S n 
for all v G A. This gives the Lipschitz property for L n as required for 
(JTB]). For ( |IcD note that certainly (x, 2n) G iS n for all x and n. Finally, 



if z G 5 n then 2; + e d G S n+1 , giving (0). □ 
Proof of Theorem [|. This is immediate from Proposition |8| above. □ 

4. Obstacles 

In this section we define obstacles, and show that they can be dom- 
inated by independent versions. Let the functions L n be defined as in 
(Q). As mentioned earlier, we say that 

site {x,2n + 1) G Z^ 1 x Z odd is good if L n (x) — 2n + 1, 

and otherwise it is bad. For y G Z^ 1 x Z CV en we define the obstacle 
at y to be 

(7) Aj, := G Z^ 1 x Z odd : y 2— »- 2 

Note that is defined only for 7/ of even height, while it consists of a 
set of sites of odd heights. 

Lemma 9. We have 

|x G Z d_1 x Z odd : x is bad j = (J 



STOCHASTIC DOMINATION AND COMB PERCOLATION 



11 



Proof. From (|5|),(|) it follows that z = (x, 2n + 1) G Z d 1 x Z odd is bad 
if and only if z G S n , which in turn is equivalent to the existence of 
y G Z ^ 1 x Z even such that y ^ yd l ^ 2 y z _ rj 

Now let (A y ) yeZ d-i xZeven be mutually independent random sets, with 
A y equal in law to A y for each y. 

Proposition 10. With the above definitions, [j A y is stochastically 
dominated by [j y A y . 

Proof. We rephrase the definition of obstacles in terms of a first-passage 
percolation model. Let each upward directed edge (z, z + e^), z G Z d 
have passage time 1 if z + is open, and if z + e<i is closed. Each 
downward directed edge (z, z + v), v G A has passage time 0. All 
other edges have passage time oo. Note that all the passage times 
are independent. We assign source time t(y) = yd/2 to each site y G 
Z d_1 x Z even , and source time oo to all sites in Z d_1 x Z odd . From the 
definition of r-open paths, for z G Z"' -1 x Z odd , we have 

(8) ze\jA y ^ T(z)<^. 

y 

Now consider a countable family of models indexed by y G Z d_1 x 
Z even . All models have the same distribution of passage times as de- 
scribed above, and are independent of each other, but in model y, the 
only source is y, with t(y) = yd/2 (all other sites have source time 
oo). Write T y (z) for the passage time to z in model y. The set of sites 
z G U 1 ^ 1 x Z odd with T y (z) < [z d — l)/2 has the same law as A y \ let us 
define it to be A y . Thus the family (Ay) has precisely the distribution 
required. Writing T(z) := mf y T y (z) and using (|), we have 

(9) ze\JA y ^ T(z)<^. 

y 

Theorem § tells us that (T(z)) is stochastically dominated by (T(z)). 
Using (^) and (^|), this implies that [j y A y is stochastically dominated 

by \J y A y , as required. □ 

5. Radii of obstacles 

Let R y be the radius of the obstacle at y, by which we mean the 
smallest r such that A y C B(y,r) (recall that B(y,r) := {z G Z d : 
1 1 2/ — 2 1 1 oo < r })- So -R y = if and only if i? y is empty. Since all 
sites in A y must have d coordinate strictly greater than y d , we observe 
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that R y is never equal to 1, and also that A y C B(y + e d , R y ). Recall 
that our geometric random variables are supported on the non-negative 
integers. 

Lemma 11. If p is sufficiently close to 1 then for each y G Z d_1 x 
Zgvenj the radius R y of the obstacle at y is stochastically dominated by 
a Geom(c) random variable, where c = c(p) — > as p — > 1 . 

Proof. Since as observed above, R y never takes the value 1, it will 
be enough to show that W(Ry > r) < c r+1 for all r > 1. We use a 
path-counting argument similar to that already used in the proof of 
Proposition 

Suppose that R y > r. Then by the definition of A y there exists z with 

||y — ^||oo > r an d y ^ d Vd — — > z. Consider some {z d — yd — l)/2-open 
A-path from y to such a z, and as before let it have U up-steps ending 
in open sites, C up-steps ending in closed sites, and D (diagonal- or) 
down-steps. Since the path is (zd — yd — l)/2-open we have 



and so C- U — D — 1 > 0. Since \\y — z^ > r, either Zd — yd > r, in 
which case U + C — D > r, or else z and y differ by at least r in some 
other coordinate, in which case D > r (since only down-steps permit 
horizontal movement). Using the earlier inequality, in either case we 
have U + C - r > 0. 

As before, let K := 3 d " 1 + 2 and q := 1—p. Then, bounding via the 
expected number of paths, 



U<(z d -y d - l)/2 — (U + C — D — l)/2, 



F(Ry > r) < 



U+C+DC 



U,C,D>0: 
C-U-D-1>0 
U+C-r>0 




X,Y,D>0 



x>o y>o d>o 



where in the second inequality we wrote X := C — U — D — 1 and 
F := U + C — r and dropped the conditions [/, C > 0. The last 
expression equals 
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(where A = A(K, q) and a = a(K, q) are defined by the last equal- 
ity), provided a < 1. Finally, we have Aa r < max(A, a) r+1 , and 
m&x(A, a) — > as q — > 0. □ 

6. Completing the embedding 

In this section we conclude the proof of Theorem |I] by combining 
the various ingredients together with some geometric arguments. We 
start by proving Lemma |7|, which states that balls may be replaced 
with sticks for the purposes of finding a Lipschitz function that avoids 
them. 

Proof of Lemma For z £ Z d we write z := (zi, . . . , so z = 

(z, Zd)- Let r > 1 (otherwise the ball and stick in the lemma are both 
empty). Suppose that {(x, h{x)) : x £ Z d_1 } does intersect B(y,r), 
say at the site u = (u,h(u)) £ B(y,r). Thus ||w — yW^ < r — 1 and 
\h(u) — yd\ < r — 1. By the Lipschitz property, the former implies 
\h(u) — h(y)\ < r — 1. Therefore \h(y) — yd\ < 2r — 2. Thus (y, h(y)) £ 
S(y, 2r — 1), and {(x, /i(x)) : x £ Z ' -1 } intersects S(y, 2r — 1). □ 

Next we check that the Lipschitz surfaces of Theorems [| and || can 
be combined to give an embedding of the comb. A slightly subtle point 
in dimensions d > 3 is that the backbone surface H will not typically 
"line up" with the stacked surfaces L n with respect to the intermediate 
coordinates 2, . . . , d — 1. Nevertheless, the use of the oo-norm in the 
definitions of ZfL gives enough wiggle room to permit an embedding. 

Suppose we are given the functions L n and H . For z £ Z d , define 
x(z) £ Z d_1 and f(z) £ 7L d by 

(10) x(z) := (z! + H((z 2 ,z 3 , . . .,z d -!,2z d + l)),z 2 , . . .,Zd-ij] 

(11) f(z):=(x(z),L 2zd+1 (x(z))). 

Lemma 12. Suppose that the functions L n and H satisfy the condi- 
tions of Theorems [| and Then the function f defined above is an 
embedding of the comb K in the open set of Z^ . 

Proof of Lemma [71]. First observe that the site f(z) is open for all z\ 
this is immediate from ([11]) and property (|Ta|) of L n . 

We next check that / is injective. From ([TdD , the sites (x,L n (x)) 
and (x',L n i(x')) are distinct whenever i / i' or n / n'. Suppose 
z ^ z'. If ,2 and 2/ differ in any of the first d — 1 coordinates then 
by (pjj|), x(z) 7^ x^ 2 )? while if they differ in the d coordinate then 
2z d + 1 ^ 24 + 1. Thus (0) gives f(z) ^ f(z'), as required. 
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To show that / is an embedding of the comb it remains to check that 

(12) ||/(z)-/(^ + e i )||oo<2 

for all z G Z d and % < d — 1, and also for i = d whenever z\ = 0. We 
first note the following key point. If Z\ — 0, then 

(13) f(z) = (x(z),2z d + l). 

This is because, by (pa|) , for w = (z 2 , £3, • • • , Zd-i, 2z d + 1), the site 
(H(u),u) is good, which means that 

L 2zd+ i{x{z)) = L 2Zd+1 (H(u), z 2 , z 3 ,..., Zd_i) 
= 2z d + 1, 



so that the two expressions in (|TTJ) and (TT3]) are the same. 

We now verify that (O) holds in the cases claimed. First suppose 
that z and z' differ by 1 in the ith coordinate, where % < d — 1, and that 
all the other coordinates agree. By the Lipschitz property ( |2"b"D of H, 
the first coordinates of x(z) and x(z') differ by at most 1, and clearly the 
same is true of the other coordinates. Hence by property (|Tb[ ), we have 
\L 2zd+1 (x(z)) - L 2zd+1 (x(z'))\ < 1. It follows that \\f(z) - /(z')||oo < 
1 < 2. 

Now suppose that z and z' differ by 1 in the last coordinate, and 
all the other coordinates agree, and suppose in addition that Z\ = 0. 
Thus 2z d + 1 and 2z' d + 1 differ by 2, and by ([13]), f(z) = (x(z), 2z d + 1) 
and f(z') = (x(z'),2z' d + 1). By (pbl) and (flO|), the first coordinates 



of x(z) and x(2') differ by at most 1, and the other coordinates agree. 
Therefore \\f(z) — /(z')lloo = 2 as required. □ 

We are ready to prove Theorem ^| and deduce Theorem |I[ 

Proof of Theorem 0. We need to show that if p is sufficiently close to 
1 there exists H satisfying fl2"a|) and (Pb]). 

Condition (^a|) says that the surface {(H(u),u)} must avoid every 
obstacle A y , and by Proposition [10], for this it suffices to instead find 
a surface avoiding the independent obstacles A y . As remarked in Sec- 
tion I we have A y C B(y + e d , R y ) (and y + e d G Z^ 1 x Z odd ), and 
by Lemma n, R y is dominated by a geometric random variable whose 



parameter c = c(p) can be made as small as desired by taking p large 
enough. Therefore it remains to show that for c sufficiently small there 
exists H satisfying (|2"b"| ) such that {(H(u),u) : u G Z d ~ 2 x Z Q dd} avoids 
lJ ygZ(i _i xZ Ad B{y,G y ), where (G y ) are i.i.d. Geom(c). Note that now 
all the relevant sites have odd heights. 
Now we map 

Z^xZodd to Z d via the transformation (m, t> , 2n+l) h-> 
(f , n, m), for u G Z d_2 and m, n G Z. It thus suffices to find a function 



STOCHASTIC DOMINATION AND COMB PERCOLATION 



15 



h : Z d_1 — > Z satisfying the same 1-Lipschitz condition (|BD as Lq, 
and whose graph {(x, h(x)) : x G Z^ 1 } avoids {J yeZ d B(y, G y ) for (G y ) 
i.i.d. Geom(c). (The transformation does not increase oo-norms). 

Now we apply Lemma [7|. The graph of h will avoid the balls B(y, G y ) 
provided it avoids the sticks S(y, (2G y — 1)+). Observe also that if G is 
Geom(c) then (2G — 1) + is dominated by a Geom(c') random variable, 
where d = y/c, so it suffices to avoid S := [J yeZ d S(y,G' y ) where (G' y ) 
are i.i.d. Geom(c'). 

The random set S consists of independent components in each of 
the lines {x} x Z, for x G Z d_1 . Within any such line, Theorem ^ 
shows that it is stochastically dominated by the open set of an i.i.d. 
percolation process with parameter c" = A^/c'. Thus the whole set 
S is dominated by the open set of an i.i.d. percolation process with 
parameter c" on Z d . Hence it follows from Theorem [5] (exchanging the 
roles of open and closed sites) that there exists a function h satisfying 
our requirements if c" is sufficiently small. Since c" = 4c(p) 1//4 , this 
holds provided p is sufficiently close to 1. □ 

Proof of Theorem [1. This is immediate from Lemma |T2| and Theo- 
rems ^] and □ 

Proof of Corollary [|. Fix k > 1 and call the site x G 7L d occupied if 
the cube kx+ [0, k) d contains some open site in the percolation model. 
For any graph G, if there exists an embedding of G in the occupied 
sites of Z^ m j then there exists an embedding of G in the open sites 

of Z| i fcm+fc _ 1 ]: we simply choose one open site from the cube of each 
occupied site in the image. Therefore, 

1 - Pc (G, Zf km+k _ 1} )\ " > 1 - Pc (G, Zf m] ). 

Setting m = 2 and G = K d , and using the fact that p c (G,Zf M i) is 
decreasing in M, the result follows from Theorem [l]. □ 

7. First-passage percolation domination 

In this section we prove Theorem (|. Recall that we have a collection 
of models indexed by i, and an additional model whose source times 
are given by inffma of source times of the others. Write Wi(e) and 
W(e) for the passage time of edge e in model i and in the additional 
model respectively. 

Proof of Theorem 0. The argument is most straightforward in the case 
where the vertex set V and the index set X are finite, and where a.s. 
the occupation times Tj(x) are finite and distinct for all i and x. (This 
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property holds, for example, when all the source times are distinct 
and finite, and each edge passage time is either oo or some positive 
continuous random variable). We begin with this case, and then extend 
to the general case by a limiting argument. 

We will define the collection of passage times (W(e)) as a function 
of the collection (W;(e)), in such a way that W(e) shares the common 
distribution of the Wj(e), that the passage times W(e) are independent 
for different e, and that T(x) < T(x) for all x. This explicit coupling 
implies the stochastic domination required. For a directed edge (x, y) 
we set 

W(x,y) := Wi(x,y), where i minimizes Tj(x). 
First we aim to show that T(x) < T(x). We have 

T(x) = mm \t(y ) + W(y k -i, y k )\. 

y m =x < ^=1 ' 

If y , yi, . . . , y m is a minimizing path in the above expression, then for 
all r with < r < m, 

r 

T(y r ) = t(y ) +J2 w (Vk-i,yk), 

k=l 

and in particular T(y r ) = T(y r _ 1 ) + W(y r -i,y r ) for 1 < r < m. 

We will show by induction that T(y r ) < T(y r ) for all such r. For the 
r = case, we have 

T(yo) = t(yo) = mmU(y ) > mmTi(y ) = T(y ). 

Now suppose T(y r _i) < T(y r _i). Let % minimize Tj(y r _i), so that 
W{y r -i,y r ) = Wi(y r -i,y r ), and T(y r ) > T(y r ) = Ti(y r ). Then 

T(y r )=T(y r - 1 ) + W(y r - 1 ,y r ) 

> 7i(y r _i) + Wi(y r -i,y r ) 

> T, t (y r ) > f(y r ), 

completing the induction. 

It remains to show that the W(e) as defined are indeed independent 
with the required distributions. We will do this by giving a different 
construction of all the models. The idea is to run them simultaneously 
in real time, revealing the random passage times only when they are 
needed. 

First consider a single model i. We begin by choosing all the passage 
times Wi(e), with the correct distributions, but we do not yet reveal 
them. (We can think of them as written on cards associated with the 
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edges, which will be turned over at the appropriate times). Label each 
vertex with a time by which it needs to be examined; initially these 
are just the source times. Now we repeatedly do the following. Find 
the vertex x with the earliest (smallest) label among those that have 
not yet been examined. Then examine x, which is to say reveal the 
passage times Wi(x, y), y G V of all edges leading out of x, and relabel 
each vertex y ^ x with the minimum of: its current label, and the 
label at x plus Wi(x,y). Repeat until all vertices have been examined. 
It is clear that the vertices are examined in order of their occupation 
times Ti(x), and that when a vertex is examined it is labeled with its 
occupation time (and this label does not subsequently change). Our 
assumptions guarantee that these times are all distinct, and so the 
choice of which vertex to examine next is always unambiguous. (These 
claims may be checked formally by induction over the vertices in order 
of their occupation times). 

Now consider simultaneously running all the models i G X in the way 
just described. We first choose all the passage times independently, 
without revealing them. At each step we examine the unexamined 
vertex with the earliest label across all the models (and we examine it 
only in the minimizing model). Clearly each individual model evolves 
exactly as before (but with its steps interspersed with the others). Our 
assumptions guarantee that no two steps are simultaneous. Finally we 
construct the passage times W(e) of the additional model: at each step, 
if the vertex that is examined (say vertex x in model i) is the first to 
be examined among the copies of that vertex x in all the models, then 
we in addition set W(x, y) = Wi(x,y) for all y G V. Since the label of 
vertex x in model % at this step is min^ T^x) (= T(x)), this agrees with 
the earlier definition of W(x, y). The key point is that the decision 
to assign Wi(x,y) to W(x,y) is made before the value of Wi(x,y) is 
revealed. It follows that the passage times (W(e)) assigned to the 
additional model are independent and have the correct distributions, 
as required. 

To extend to the general case, we consider a sequence of approximat- 
ing finite systems of the kind just considered. Without loss of general- 
ity, suppose that the index set X is contained in N. In the nth system in 
our sequence of approximations, we take a finite vertex set V^ n \ such 
that V {n) | V as n — > oo. All source times and passage times involving 
a vertex not in are set to infinity. Any source time for a vertex in 
that was previously be set to infinity is now instead given value n 
(to ensure that every vertex in is reached in finite time). Further- 
more we consider only models indexed by i G {1, 2, . . . , n}. Finally we 



18 



ALEXANDER E. HOLROYD AND JAMES MARTIN 



perturb the source times of vertices in V^ n \ and the passage times of 
edges joining points of V^ n \ by adding an independent Uniform(0, -) 
random variable to each. (This means that the source times are no 
longer deterministic - however, we can regard the randomness as being 
on two levels: given any choices of the source times, we have a set of 
models with random passage times). This ensures that a.s., the finite 
system satisfies all of our earlier assumptions. 

Write T[ n \x), f^ n \x) and T^(x) for the passage times in the nth 
approximation. These quantities are finite for any x G V n . But also, 
since any set of vertices is eventually contained within V n , and each 
model i is eventually included in the system, it follows that for any 
finite set A, the random vectors (T^ n \x)) xG A and (T^ n \x)) xe A converge 
in distribution as n — )■ oo to (T(x)) xe A and (T(x)) x& a respectively. We 
know from the argument applied to the finite case that (T^ n \x)) xe A 
is stochastically dominated by (T^ n \x)) x€ A for any n. Hence from 
the convergence in distribution as n — > oo, we obtain also that in 
fact (T(x)) x( za is stochastically dominated by (T(x)) xe A- Since this 
holds for any finite subset A C V, it follows that in fact (T(x)) x& y is 
stochastically dominated by (T(x)) xE y, as desired. □ 

We remark that Theorem |] may easily be extended for example to 
models with undirected edges instead of (or in addition to) directed 
edges, or with passage times at sites. As long as all the passage times 
are independent, exactly the same methods used above will continue 
to apply. 

8. Domination in one dimension 

In this section we prove Theorem |3|. We start with the following 
one-sided version. The interval [a, b) is taken to be empty if a = b. 

Proposition 13. Fore G (0,1), let (G n ) rig z be i.i.d. Geom(c) random 
variables. The random set Z R U n .ez[ n ' n + G n ) is stochastically dom- 
inated by the open set of i.i.d. site percolation on Z with parameter 
min(2c, 1). 

Proof. Define the indicator variable 

Bi :— 1 \i G [n, n + G n ) for some fiGZ] 

= 1 \G n > i — n for some n < i \ . 

Then we must prove that [Bi)i & z is dominated by an i.i.d. Bernoulli(2c) 
sequence (when 2c < 1). For this, it is enough to show that a.s., 

(14) ¥(Bi = 1 | (Bj)^ < 2c 
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for all i (see e.g. [[TJ, Lemma 1]). Since the process (Bi) is stationary 
it suffices to show ( |H|) for i — 0. 

We may think of the system as an M/M/oo queue in discrete time. 
At each time n, a customer arrives whose service time is G n . The 
customer will depart at time n + G n , and so occupies the system during 
the interval [n, n + G n ). (If G n = then the customer is never seen at 
all). Now Bi is the indicator of the event that there is some customer 
present at time i. 

We introduce the key random variable 

N := max{n > : G_ n > n}. 

We can think of N as the age of the oldest customer who has not left 
the system before time 0. (For example, if iV = then all the previous 
customers have left before time 0; on this event we have Bq = 1[Gq > 
0], since the only customer who could be present at time is the one 
arriving at time 0.) The Borel-Cantelli lemma shows that iV is a.s. 
finite. 

We claim that a.s. 

(15) P(B = 1 | (B 3 ), <0 ,N) = F(B = l\N). 
To prove this, we must establish 

(16) F(B = 1 | E, N = n) = F(B = 1 | iV = n), 

for all events E G cr((i?j)j<o) for which the conditional probability on 
the left exists. Observe first that N = n forces _B_ n = ■ ■ ■ = B-i = 1, 
so it is enough to prove flTBD for E G a^B^j^n). 

To verify the above, observe that the two families C := {Gj)j < - n 
and 1Z := (Gj)j>_ n are conditionally independent of each other given 
N = n: this is because they are independent without the conditioning, 
while {iV = n} is the intersection of an event in cr(£) and an event in 
a(1Z). Now, [Bj)j < _ n is a function of C. On the other hand, on N = n, 
we have that Bq is a function of 1Z (since N = n guarantees that any 
customer arriving before time — n has already left the system before 
time 0). It follows that (Bj)j < ^ n and Bq are conditionally independent 
given N = n; this gives precisely fllTf ) for E G cr((5 J -)j < _ n ) as required, 
thus proving the claim (|TJ). 
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Returning to the proof of (14), we have 
¥>(B = 1 | JV = n) 

= W(G-j > j for some j > | G_„ > n, and G_j < j for all j > n) 
= P(G_j > j for some < j < n | G_ n > n) 
= P(G_„ > n | C_ n > n) 

+ [1 - P(G_ n > n | G_ n > n)] P(G_ j > j for some < j < n) 
< c + (1 - c)(c + c 2 + • ■ ■ + c n ) < 2c. 

Combining with flT5|), we have shown that F(Bq = 1 | (Sj)j<o, N) < 2c 
a.s.; then averaging over N gives (0) (for i = 0), as required. □ 

Proof of Theorem^. Let (Ajnez an d (-Rn)ngz be i.i.d. Geom( v /c) ran- 
dom variables. Let G n := min(L n , R n ), which is Geom(c). Then, with 
all intervals understood to denote their intersections with Z, 

\J(n-G n ,n + G n ) C (|J (n - L n ,n\j U ^|J[n,n + i? n )J. 



By Proposition [13| and symmetry, the right side is dominated by the 
union of the open sets of two independent percolation processes each 
of parameter 2^c, which is itself the open set of a percolation process 
of parameter 1 — (1 — 2^/c) 2 < 4-^/c. □ 
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